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In this paper we study analytically a simple one dimensional model of mass transport. We 
introduce a parameter p that interpolates between continuous time dynamics (p ^ limit) and 
discrete parallel update dynamics {p — 1). For each p, we study the model with (i) both continuous 
and discrete masses and (ii) both symmetric and asymmetric transport of masses. In the asymmetric 
continuous mass model, the two limits p — 1 and p —> reduce respectively to the g-model of force 
fluctuations in bead packs [S.N. Coppersmith et. al., Phys. Rev. E. 53, 4 673 (1996)] and th e 
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^^ ■ recently studied asymmetric random average process [J. Krug and J. Garcia, 3ond-mat/9909034| 

fT^ ' We calculate the steady state mass distribution function P{m) assuming product measure and show 

that it has an algebraic tail for small m, P{m) ^ m~^ where the exponent /3 depends continuously 

_i on p. For the asymmetric case we find /3(p) = (1 — p)/{'2 — p) for < p < 1 and /3(1) = —1 and for 

"y^ \ the symmetric case, /3(p) = (2 — p)^/(8 — 5p + p^) for all < p < 1. We discuss the conditions under 

flj ■ which the product measure ansatz is exact. We also calculate exactly the steady state mass-mass 

r^ ' correlation function and show that while it decouples in the asymmetric model, in the symmetric 

" case it has a nontrivial spatial oscillation with an amplitude decaying exponentially with distance. 
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I. INTRODUCTION 



There is a wide variety of physical systems in nature where the basic microscopic dynamical processes involved 
■ O ■ are aggregation, fragmentation, adsorption, desorption and transport of mass. These processes are abundant and 
occur in systems such as colloidal suspensions [nl, polymer gels [pl0], river networks B, aerosols and clouds [p| and 
surface growth phenomena involving island formation p] . These systems can have different types of non-equilibrium 
stationary states and phase transitions between them as the rates of the underlying microscopic processes are varied. 
While for systems in thermal equilibrium the stationary state is characterized by the Gibbs measure, there is no such 
general recipe for non-equilibrium systems. In order to gain more insights on the nature of these steady states and 
\^ • possible phase transitions between them, several simple lattice models involving mass transport have been proposed 
^^ ' and studied recently mM ■ By virtue of their simplicity, these lattice models are often amenable to exact analysis and 
Cn ■ yet contain rich and nontrivial physics. 

These models constitute simple examples of interacting many particle systems out of equilibrium; in particular the 
dynamics of these systems do not obey detailed balance. The steady states of interacting many body systems are in 
^s^ ' general difficult to characterize and only a few exact results are available. These include simple exclusion processes 
"Y*^ with open and closed boundary conditions |10|] , abelian sandpile models of self organized criticality |ll| , trafhc models 
Cd [Q and mass aggregation model of Takayasu . Moreover the steady states in some cases are non universal and 

H depend on the detailed nature of the dynamics used for updating. For example the steady state in the asymmetric 

' [ simple exclusion process depends on whether the update rules are parallel or random sequential [[l3| . It is therefore 
^ • desirable to study more of such simple models in a systematic way in order to get insight into the nature of the 
Q non-equilibrium steady states. In this paper we study a simple lattice model of mass transport analytically which 
O ' sheds some light on these general issues pertaining to interacting many body systems. 

^ I Besides, the dynamics in seemingly unrelated systems can often be mapped onto simple one dimensional mass 

• "^ . transport models evolving with time according to some prescribed rules. These systems include river networks 

rN iQjQ 1 force fluctuations in granular systems such as bead packs jl^] , traffic flows |1^ , voting systems jl^ill^ > wealth 

S . distributions ||l^, generalized Hammersley process iQ a-nd inelastic collisions in granular gases pi[ |. 

■ ■ ■ ' In this paper we study analytically a model of mass transport in a one dimensional lattice. Each lattice site contains 

a nonnegative mass variable and the dynamics consists of transporting a flnite amount of mass from each site to its 

neighbours. The amount to be transported is randomly chosen from a given distribution. We introduce a parameter 

p that interpolates between continuous time dynamics (p — *■ limit) and discrete parallel update dynamics {p — 1). 

For each p, we study the model with (i) both continuous and discrete masses and (ii) both symmetric and asymmetric 

transport of masses. 

The paper is organized as follows. In section II, we deflne the mass model precisely and discuss its mapping to 
other models of non equilibrium statistical mechanics. We also summarize our main results. In section III, we discuss 



the asymmetric continuous mass model and solve for the mass distribution function in the steady state assuming that 
product measure holds. We then discuss under what conditions the product measure is exact. The two point mass 
correlation function is also computed exactly. In section IV, we study the symmetric version of the continuous mass 
model. We show that product measure becomes exact in a particular limit. We also compute the stationary two point 
mass correlation function exactly for the symmetric model. In section V, we study the discrete mass version of the 
model. Finally we conclude with a summary and outlook in section VI. Appendix A contains a proof that product 
measure fails for the asymmetric model for any p < I. In appendix B we prove that product measure is exact for the 
symmetric model in the p — > limit. 

II. THE MODEL 

Our mass model is defined on a lattice. For simplicity, we define it here on a one dimensional lattice with periodic 
boundary conditions. The generalization to higher dimensions is straightforward. Each lattice site contains a nonneg- 
ative mass variable. We consider two versions of the model: (i) when the mass at each site is a continuous variable 
and (ii) when the mass at each site is discrete. 

First consider the continuous mass model (see figure 1). We start with a random configuration of masses, m{i) at 
each site i. The dynamics is defined as follows. Each discrete time step consists of two moves: (1) a fraction iji of each 
mass rrii is chosen and then (2) with probability p, this fraction is added to a neighbouring site and with probability 
(1 — p) it remains at the original site. In the asymmetric model, the fraction is always added to the right neighbour. 
In the symmetric case, the fraction goes either to the left or the right neighbour with equal probability. We note that 
total mass is conserved by the dynamics. Thus the model has two parameters, the probability p and the average mass 
per site p. 

The fractions rji are independent and identically distributed random variables drawn from a probability distribution 
on [0, 1]. In this paper, we will mostly consider a uniform distribution of r/i on [0, 1] though some of our results can 
be generalized to a class of other distributions. 

We note that for p — 1, all the chosen fractions of masses are definitely transported to their neighbours. This 
corresponds to fully parallel update dynamics. In this case, the asymmetric version of the model reduces exactly to 
the g-model introduced by Coppersmith et. al. (l5| to study force fluctuations in random bead packs. In this case, 
the mass at each site evolves as 

mi{t + 1) = 77i_imi_i(t) + (1 - rii)mi{t). (1) 

In the context of bead packs, the indices i and t index the site i at a depth i in a two dimensional packing. Then 
for large t, rrii represents the force supported by a bead at (j,i) scaled by the mean weight and rjirrii is the random 
component of the weight (scaled by the mean weight) transmitted from a bead at depth i to a neighbouring one at 
depth (i + 1) that touches it. The same equation was also studied in Ref. M in the context of a lattice gas model. 
The stationary mass distribution P[m) of the q-model was solved exactly ||l5|| and remarkably the mean field theory 
turned out to be exact in the thermodynamic limit for the case of uniform distribution of the fractions 77.; 's. This 
means that the steady state joint distribution of masses at different sites factorises, P(r7T,i, 7712,7713 . . .) = Ili^l"^')- 
In other words, the product measure is exact in this case and P{m) was shown |L5|,p[ to have a simple distribution, 

P{m) = ^6-2™/". (2) 

P 

In the opposite limit p — > 0, the probability that two or more sites will be simultaneously updated in a single 
move is 0{p^) and hence negligible. With the choice of p = Ai, this limit thus corresponds to the random sequential 
continuous time dynamics. This case has been studied recently by Krug and Garcia [1^ and assuming that product 
measure holds they derived the steady state mass distribution P{m) for uniform distribution of 77; 's, 

P{m)^ —=^=ey.Y>{-m/2p). (3) 

^ iTipra 

The difference in the small 777 behaviour of the mean field Pirn) in the parallel and random sequential case was 
correctly noted by Krug and Garcia pq | . We have studied both the symmetric and asymmetric versions of the model 
for arbitrary p. Our main results are summarized as follows: 

(1) In the asymmetric case, we show that within the mean field approximation P{m) ^ m~^ for small tti, while 
it decays exponentially for large m. The exponent I3{p) depends continuously on p with a discontinuity at p = 1, 



P{p) = (1 — p)/{2 — p) for p < 1 and /3(1) = —1. In the symmetric case, the mean field P(rn) has a similar behaviour 
except the exponent /?(p) = (2 — p)^/(8 — 5p + p^) for all p in [0, 1]. Note that in the symmetric case, there is no 
discontinuity dX p = 1 . 

(2) In the asymmetric case, we prove rigorously that the product measure is exact only for p — 1. For any 
p < 1 (including the random sequential p = case), we show that the product measure ansatz, P(mi, m2, m.3, ■ • •) = 
Y[i Pirrii) breaks down. But remarkably the mean field P{m) is almost indistinguishable from the P{m) obtained from 
numerical simulation in one dimension. We note that the breakdown of product measure property does not necessarily 
mean that the correct single point distribution P{ni) is still not given by the mean field P{m); in fact numerical results 
strongly suggest that the mean field P{m) is exact even though product measure fails. In the symmetric case on the 
other hand, product measure is exact only for p — > but fails for any p > 0. Besides, as opposed to the asymmetric 
case, the mean field P{m) is considerably different from the distribution obtained numerically. This is due to strong 
correlations between masses in the symmetric case as mentioned below. 

(3) The two-point mass correlation function C{r) = {m(0)m{r)) between two sites at distance r can be computed 
exactly (without recourse to the assumption of product measure) for arbitrary p in both asymmetric and symmetric 
models. We find that in the asymmetric case, for all values of p, the connected part of the correlation function 
vanishes, C{r) — p^ = for r > 0. In the limit p ^ 0, this fact was noted by Krug and Garcia ||l6[. This however does 
not imply the validity of product measure is exact which would require factorization of all higher order correlations 
as well. For the symmetric case, the correlation function factorises only for p = 0. However for p > 0, the function 
C{r) — p^ has a nontrivial spatial dependence. It oscillates with distance r and the amplitude of the oscillation decays 
exponentially with r. 

We have also studied a discrete mass version of the above model. In this case the mass nii at any site i can take 
only discrete non-negative integer values. Instead of a random fraction breaking off a mass as in the continuous case, 
the mass to be taken out of a site is a random variable that takes only discrete values 0,1,2,3 . . . , rrii equally likely, 
i.e., any of these values is chosen with the same probability l/{mi + 1). Then as before, with probability p, the chosen 
mass is actually transported to a neighbouring site and with probability (1 — p) it stays at its original site. We derive 
the explicit expressions for the mass distribution for the discrete case also. 

This model can be mapped onto a model of hard core particles moving with long range jumps in one dimension 
MEm. First consider the continuous mass model. Each site of the lattice corresponds to a particle (point) on the real 
line and the mass rrii represents the continuous gap between i-th and (i -\- l)-th particle. The transport of random 
fraction of m^ from the «-th site to [i + l)-th site corresponds to the {i + l)-th particle jumping to the left by a random 
fraction of the available gap between it and its left neighbour (see figure 2). The discrete mass problem similarly 
corresponds to particles moving on a one dimensional lattice (as opposed to the real line in the continuous case) with 
hard core repulsion. At each time step a particle moves to a site randomly chosen from the set of empty sites in front 
of it. This is a generalization of the simple exclusion process where a particle can jump only to a nearest neighbour 
site provided it is unoccupied. In this generalized case, while the hard core repulsion is respected, long range jumps 
are allowed. 

The discrete mass problem can also be mapped onto an interface growth problem via the usual mapping from a 
lattice gas model to a growing interface. Starting from a reference height h — 0, a particle at site i corresponds to 
h{i + 1) = h{i) — 1 while a hole corresponds to h{i + 1) = h{i) + 1. Under this mapping, our problem corresponds to 
the following rules: Any stretch of the interface with slope equal to 1 can be split at any randomly chosen point in 
between into two sections of slope 1 connected by a bond of slope —1. 

III. THE ASYMMETRIC MODEL WITH CONTINUOUS MASS 

In this section, we study the continuous mass model where in each time step a fraction of the mass from any given 
site is transported with probability p to its right neighbour. The Langevin equation for the mass update can be 
written as 

■mi{t + 1) = mi{t) - airiimi{t) + (7i^i?7i_imi_i(i), (4) 

where the fractions ?7i's are random numbers in [0, 1] chosen from a uniform distribution and the random variables 
ai 's take values 1 with probability p or with probability q = 1 — p. The distribution of both of these variables are 
independent from site to site. Defining r^ — rjiUi, we get 

mi{t + 1) = mi{t){l - Ti) + mi_i(t)ri_i, (5) 



and it is not difRcult to see that the effective distribution /(r^) of the random variable r^ on [0, 1] given by 

/(rO = g<5(r,)+p. (6) 

The evohition equation of the single point mass distribution function P(mi,t) (which is independent of i due to 
translational invariance) can be written down exactly, 

P{mi,t+1) ^ / drrii^i / (irj_i/(rj_i) / dm^ / drif{ri)P{mi^i,m[,t) 
Jo Jo Jo Jo 

X S{m'i{l- ri) +mi^iri-i - nii). (7) 

Here P(?7ii_i , m^, t) is the joint probability distribution of mass rrii-i at site i — 1 and rrii at site i. The time evolution 
of the single point probability distribution involves the two point joint probability distribution function. Similarly the 
n-point probability distribution will involve the (n+l)-point joint probability distribution and in general this hierarchy 
cannot be broken. 

A. Mean Field Theory 

We first compute the single point mass distribution P{m) from Eq. (|7[) by assuming that the joint distribution 
factorises in the steady state, P(?7ii_i, rrii) = P{jni-i)P(mi). This approximation clearly ignores correlations between 
masses. Within this approximation, Eq. (|7|) involves only single point distribution function P{m). Taking the 
stationary limit, t — > oo, and using the explicit form of the distribution /(r) from Eq.W), we find that the Laplace 
transform, Q{s) — L P{m) exp{~ms)dm satisfies the equation, 

where V{s) — L Q{sr)dr. We note that J^(s^) = Q{s)- Eliminating Q, we get a first order differential equation for 
V which can be integrated to give, 

^ " ^ = W2, (9) 



y2-p 



where the integration constant has been determined by using the fact that dQ / ds\g^Q — —p with p being the average 
mass per site. The above equation reduces to a quadratic, cubic and linear equation in V for p — Q, p — 0.5 and p = 1 
respectively. 

For the fully parallel dynamics p = 1, we get V{s) — 2/(2 + ps) and hence Q{s) — 4/(2 + ps)^ . By inverting the 
Laplace transform, we recover the result P{m) — inig^^™''' obtained by Coppersmith et. al |lj]. In the random 
sequential limit, p -^ 0, we get 



V{s) = 


ps 


Q{s) = 




x/1 + 2ps 


°(™) = 


_ exp( m/(2p)) 



(10) 

y/ Z,7l pill 

The same result was obtained by Krug and Garcia n& by a somewhat indirect method by computing the moments 
and then guessing the distribution from its moments. When p = 0.5, the expression for P{m) is complicated and we 
do not reproduce it here. 

For arbitrary p, a closed form expression of P{m) is difRcult to obtain. However the asymptotic behaviour of P{m) 
for large and small m can be easily derived. For large to, we expect P(rn) '^ e~"™. The decay coefficient a can be 
derived by noting that the Laplace transform Q{s) must have a pole at s = —a. From Eq. m), we note that the pole 
of Q occurs when V = {1 — q^)/pq. Using this in Eq. (^, we get, 

2{1-pY-p nn 



In the limits, p = I and p ^ 0, this gives the correct decay coefficient of P{rn). 

For small to, on the other hand, P(rn) has an algebraic tail, P{ni) ^ m~^ . From Eq. (O), we note that for large s, 
V{s) w (2/ps)i/(2-P). Using Q{s) = d{sV)/ds, we get for large s, 

Q{^)^\r^{-Y'''-'^- (12) 

2 — p ps 

This implies that for p < 1, P{m) ~ m"^ for small ?ti with /3 = (1 — p)/(2 — p). Note that for p ~ 1, the coefficient 
of 1/s vanishes in Eq. (|l3) and the leading order term decays as 1/s^, implying (3 — —1 for p = 1. Thus there is a 
discontinuity in the exponent f3{p) at p = 1. 

How good is the product measure ansatz? In general, we have noted before that the equation of the n-point 
distribution function contains the (n+l)-point distribution function. If the product measure ansatz were to be exact, 
then one has to ensure that every equation of the hierarchy is satisfied by the product measure ansatz. This was 
in fact proved to be case for p = 1 fl^]. It is easy to show that this ansatz is exact only for p = 1 and fails for all 
p < 1. This is proved by showing that for p < 1, the second equation of the hierarchy (involving the two-point and 
three-point distributions) is not satisfied by the P{m) obtained from the ffist equation of the hierarchy, i.e., Eq. (0) 
assuming product measure. Algebraic details are given in Appendix A. 

For p < 1, we compared the mean field answer for P(m) with the numerically obtained distribution in one dimension. 
In the limit p — > 0, the mean field P{in) matches extremely well with the numerically computed one. This was also 
noted by Krug and Garcia |16| . For arbitrary p, we do not have a closed form expression of mean field P{m) to 
compare with the simulation results. However, various moments of to with the mean field P{m) can be computed 
exactly for arbitrary p and compared to the numerically obtained moments. The mean field moments are computed 
by expanding V{s) in powers of s. We list the the moments < m" > for n = 1, . . . 5 below. 

< m > = p 

< to3 > = 3(2-p)^(5-3p)^3 (13) 

4 5{2-p){2l-26p + 8p^) 4 

< TO > = p 

5 _ 15(2 - p) (504 - 955p + 600p^ - 125^^^) 5 
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To check how accurate these mean field moments are, we have computed these moments directly from numerical 
simulation of the model for different values of p. In Table I, we compare the mean field moments (up to order 5) to 
the numerical ones for a representative value of p = 0.8. The closeness of these moments to the numerical values for 
arbitrary p suggests strongly that the mean field P{m) may be exact for all p even though the product measure fails 
for p < 1. 

B. Correlation Function 

In this subsection we compute the two point mass correlation function exactly for the asymmetric continuous mass 
model. In the continuous time case (p ^ limit of our model) , this was computed exactly by Krug and Garcia pq ] 
for arbitrary probability distributions of the random fraction r. We reproduce their calculation here for completeness. 
Multiplying mi{t+ 1) by mj{t+ 1) in Eq. M) and taking expectation value in the steady state, we find that two point 
correlations Cj =< miimi^j > satisfy the following set of linear equations, 

Co(Mi - M2) - Ci^ii{l - fii) = 0, 
Co(/ii - Ai2) - 2CiMi(l - /^i) + C2Mi(l - /ii) = 0, (14) 

J-2,3,4,... 

with the boundary conditions Cj -^ p^ as j — > 00. Here /ii = (r^) and p,2 = {ri^) are the first and second moments of 
the random fraction r^ distributed according to Eq. (0) . These set of of equations can be solved easily to give |lq] , 





Co(mi - 


M2) -Ciiiiil - 


-/-^i) 


= 0, 


^^2)- 


2CiMi(l- 


/^i) + C2Mi(l - 


-/^i) 


= 0, 






Cj-i — 2Cj + Cj+i 


-0, 


-p' 


as j — > 00 


. Here /ii = (r^) 


and 


M2 = 



_ Ml(l- AH) „2 
Oq — p 

Ml -M2 
Q=p2, J = 1,2,3,... (15) 

Thus the two point correlation function < miinrij > is equal to < rrii >< rrij > for i ^ j. In fact this conclusion holds 
for any arbitrary distribution of the fractions r^. This however does not mean that the product measure is exact. 
That would require that all higher order correlations must also factorize. In fact, for the asymmetric model, it can 
be shown [E2| that the 3-point correlation function does not factorize except for p — 1. 

IV. THE SYMMETRIC MODEL 

In this section, we study the continuous mass model where mass transport has no bias in direction. Once again we 
have a continuous mass rrii at each site. In each time step, a fraction is chosen at random from a uniform distribution 
on [0, 1] and this fractional mass is transported to the left or right nearest neighbour with equal probability p/2. 
With probability g = 1 — p, the fractional mass stays at the original site. In order to write down the mass evolution 
equation, it is convenient to define a set of variables Sj at each site i. The variable Si can be either +1 or —1 with 
equal probability 1/2. If Si — 1, it indicates that the fractional mass from site i is transported to the right neighbour. 
On the other hand, Si = — 1 indicates transport to the left neighbour. Then the mass evolution equation can be 
written down as in the asymmetric case, 

mi{t + 1) = (1 - ri)mi{t) -\ -^^^n^imi^i{t) H -^ — ri+im^+i, (16) 

where the random variables r.i have the same distribution /(r^) = qS{ri) +p as in the asymmetric case. The evolution 
of the single point mass distribution function P{m, t) can be written down as in the asymmetric case (Eq. (|7|)). The 
only difference is that now the single point equation contains three point distribution (as opposed to the two point 
function in the asymmetric case) and the additional Si variables. 

A. Mean Field Theory 

Assuming product measure, this equation can be solved in the same fashion as in the asymmetric case. It follows 
that the Laplace transform Q{s) of P{m) in the stationary state satisfies. 



pV{q+pV + l)^ 
i-q{q + pV+iy 



Qi^) - : ""J'lJrL (17) 



where V{s) = /„ Q{su)du as in the asymmetric case. Using Q{s) — d{sV)/ds, we find the function V{s) is given by 
the solution of the following nonlinear equation, 

[l-^^]^^-^^^(l-^) = f. (18) 

In the limit p — > (random sequential limit), this equation can be solved in closed form and we get, 

Pim) = -j^= exp(-m/(2p)), (19) 

which has the same expression as for the asymmetric p — > case. For other values of p, while we are unable to 
get a closed form expression, the asymptotic behaviour of P(m) for large and small m can be easily derived. We 
find that for large m, P{m) ~ exp(— am) where the coefficient a{p) can be determined in the same way as in the 
asymmetric case. It is given by a long expression which we do not present here. For small tti, P{m) has an algebraic 
tail, P{m) ~ m~^ where the exponent (3{p) can be determined by examining the large s behaviour of Q{s). We find 
(i{p) = (2 — p)^/(8 — 5p + p^) which decreases continuously from 1/2 [p -^ 0) to 1/4 {p = 1). 

For the symmetric case, we show in Appendix-B that the product measure is exact in the p — > limit. For p > 0, 
the product measure fails and unlike the asymmetric case, the mean field P{m) is considerably different from the 
distribution obtained numerically. This failure of mean field theory for p > Q shows up in the calculation of two point 
correlation function as done in the next subsection. However, while the mean field theory fails for large m (as evident 
from expectation value of the moments of the mass distribution shown in Table II), it matches very well with the 
numerical result for small m (see figure 3). 



B. Correlation Function 



For the symmetric model, the translationally invariant stationary two point mass correlation function, Cj-i = 
{mirrij) does not factorize for j ^ i as in the asymmetric case. Below we compute the two point correlation exactly 
and show that the connected part of the correlation function in fact has a nontrivial spatial dependence. 

Multiplying Eq. {nm by mj{t + 1) and taking expectation value, we find that in the stationary limit t — > oo, the 



correlation function Cj satisfies, 

MiQ-2(l - S,2) + 4(1 - /ii)Q-i + 2(3/ii - 4)Q + 4(1 - /ii)C,+i + /iiQ+2 - j - 2, 3, . . . 

4(1 - w)Co + 2(7//i/2 - 4)Ci + 4(1 - //i)C2 + M1C3 = (20) 

4(-l + w)Co + 4(1 - /^i)Ci + M1C2 = 0, 

where /^i = (r^) and fi2 = {fi^) are respectively the first and second moments of /(rj) and w — 112/1^1- 

Let G{z) — X]^i C'j-^'' be the generating fmiction. Multiplying Eq. ( po[ ) by z^ and summing over j's, we get 

^^'^ - (l-z)[4z + /ii(l-z)2] ^21) 

The boundary condition, Cj -^ p^ as j ^ 00 implies that G{z) -^ /5^/(l — z) as z -^ 1. This gives us one relation 
between Co and Ci, 

Ci = (2p2_2(l-w;)C'o)/m- (22) 

We need one more condition to fix both Cq and Ci. This is obtained by noting that G{z) in Eq. ( pl| ) has three poles, 
z = 1 and z = z± where z± = {fj,i — 2 ± 2-^1 — /zi)//Lti. We note that jz+j < 1 which would imply that Cj will blow 
up exponentially as \z+\^ for large j. Since this can not happen, the numerator on the right hand side of Eq. ( [2l| ) 
must also vanish at z — z^ in order to cancel the pole. This provides an additional condition which together with Eq. 
(H) gives, 

_ p'^ii + z+) _ p^yi-pi 

(1 — u;)(l — z+) 1 — w 



and Ci can be determined from Eq. (^21). Inverting the generating function, we find that for any n > 0, 

Cn^p'[l-zl]. (24) 



Since z+ = (/ii — 2 + 2^/1 — j.ii)/fii lies in [—1,0], clearly the connected part of the correlation function has a nontrivial 
oscillation with distance with an amplitude that decays exponentially with the distance. 

Curiously the function d for n > depends only on fii but not on /X2, whereas Co involves both pi and p2- We 
also note that the above exact result is valid for any arbitrary distribution /(r) of the fractions n and not just for the 
special distribution given by Eq. (|6|). For that distribution, we get from Eq. (|6|), pi — p/2 and p2 — p/3 and hence 
w — 2/3. One useful check is that in the limit p ~f 0, we get z+ — *■ implying complete decoupling of the two point 
correlation. This is consistent with the fact that product measure is exact in the symmetric case only in the p — > 
limit. For p > 0, the correlation function has a nontrivial spatial dependence and product measure clearly fails. 

V. THE DISCRETE MASS MODEL 

In this section we study the model when the mass rrii at each site i is a discrete non negative integer. In each time 
step, a block of size Ui is chosen at each site and is transported to its neighbour with probability p and stays at the 
original site with probability q = I ~ p. The block size Ui is a discrete random variable which can take values 0, 1, 
2, . . ., rrii, all with equal probability l/{mi + 1). As in the continuous mass model, the mass transport can be either 
asymmetric or symmetric. We study here only the asymmetric model but the symmetric version can be studied by 
using similar procedures. 

There is an equivalent lattice gas representation of this model in one dimension as mentioned in Section- (II). In 
this mapping, lattice site i of the mass model corresponds to the i-th hard core particle and the mass mi represents 
the number of holes or empty sides between the i-th and {i + l)-th particle. In the lattice gas dynamics of the 



asymmetric model, at each time step every particle jumps to any one of the available vacant sites in front of it with 
equal probability. 

The analysis of the stationary mass distribution of the asymmetric discrete model proceeds along the same line as 
its continuous counterpart. We write down the evolution equation of the single site distribution function P{m, t) in 
terms of the joint two point distribution P{mi,m2,t). Assuming product measure holds, the evolution equation is 
given by, 

■^-^ ^-^ ^-^ ^-^ (rrii-i + 1) m + 1) 

m._i=0mi=0m;=0m2=0 ^ '^ ^ ' 

^-^ ^-^ m' + 1 ) 
+ M E E E ^p-^^Sim'^ + m,-m.) 

mi-i—0 7ni—0 rii.—O 

oo 

+ g2 ^ F(m:)(5(m: - m,). (25) 

We define the generating function, Q{x) = J^'o' P{m)x"^. In the stationary limit, we get from the above equation, 

^, ^ _ [f{x)-f{l)Mf{x)-f{l))+q{x-l)] 

where f{x) = 'J2m=o m+i ^"^^^ ■ Using Q{x) — -J^, one can obtain closed form expressions of Q{x) and hence of P{rn) 
only in the two limits, p — 1 and p — > 0. For the fully parallel dynamics {p — 1), we get, 

«(^) - (i-f(Li))^ (^^) 

_ 4(m+l)p" 
^^ " (p + 2)"+2 ■ ^^^' 

For the random sequential case (p —^ 0), we find 

^l-2p(a;- 1) 

^^" (l + 2p)»+l/2(„!)222n- ^-^^^ 

It is easy to check that in the limit of large m and p, these distributions reduce to their continuous counterparts (Eq. 
(0) and Eq. (||) respectively) as expected. 

As in the continuous asymmetric model, it turns out that the product measure is exact only in the p = I limit. 
The proof that the product measure is exact for p = 1 in the discrete case can be derived by following the same line 
of arguments as used for the continuous case ||I^. Basically, one writes down the exact evolution equation for the 
n-point joint distribution which involves the (n + l)-th point joint distribution. One makes the ansatz for product 
measure and ensures that this ansatz is consistent for all values of n, i.e., all the equations of the hierarchy satisfy the 
product measure ansatz. 

Without giving the details we just outline below few basic steps. For p = I, assuming product measure in the 
equation involving single point and two point distributions (Eq. Eq), we obtain P(rn) as given by Eq. ES. Consider 
first a cluster of n neighbouring sites 1,2, .., n. The time development for the n- point probability distribution can be 
written as 

oo mo oo mj, tdI I I I j-\ 

P(mi,...,m„,t + 1)= E E--- E E P(^o,m„...,m„,t) 



-^ (m'r. + l)(m', + 1) . . . (m' + 1) 
(5(mi - ri + ro - mi) . . . b{m^ - r„ + r„_i - m„). (31) 



We have to now show that the product measure ansatz P(mi, TO2, ■ ■ ■) = Yii P{''^i) in the steady state with P{m) given 
by Eq. ( [28| ) is consistent with Eq. (|3^). To show this, we consider the n- variable generating function Q{xi, . . . , a;„) = 
J2m =0 ■ ■ ■ J2m =0 ^(™i' ■ ■ • I Wri)a;i'^ • ■ ■ 2;™". Wc assume product measure on the right hand side of Eq. (^), sum 
over the rrii^s and obtain, 

^f N _ a- fjXl) f{x2) - fjXl) fjXn) - f{Xn-l) a - .fjXn) 

(^[Xi, . . . , Xn) — • ■ • 

1 X± X2 X± Xn Xn—l -L X^ 

= Q{xi)Qix2)...Qixn) (32) 

where in deriving the last step we have used Q{x) — df /dx and the expression of Q{x) from Eq. (J2q). One can repeat 
the same calculation when the n sites are not necessarily neighbours. This therefore proves that every equation of the 
hierarchy of distribution functions satisfies the product measure ansatz for p = 1. However, this proof fails for p < 1 
as in the continuous case and the same line of argument used for the continuous case (see Appendix A) goes through 
for the discrete case. Even though the product measure fails, the mean field answer for other values of p match very 
well with the numerically computed one. For the random sequential case (p — > 0), we compare the mean field result 
for the single site probability distribution with numerical simulation (see figure 4). 

VI. SUMMARY AND CONCLUSION 

In this paper we have studied a simple mass model of chipping and aggregation where a mass at a site can chip off 
a fraction to its neighbour. A parameter p was introduced which allowed us to interpolate between parallel dynamics 
and random sequential dynamics. We studied the model for both continuous and discrete masses as well as for 
symmetric and asymmetric transport of mass. 

We have calculated analytically the mass distribution function P{ra) in the steady state for all p assuming product 
measure, i.e., neglecting correlation between masses. In some cases we proved that this product measure is exact. 
One of the main results is that the distribution P{m) has an algebraic tail for small to, P(m) ~ mT^'^^^ where the 
exponent (i{p) depends on p. Thus the steady state is non universal and depends on the specific nature of the dynamics 
characterized by the parameter p. 

Another interesting point is that for the asymmetric continuous mass model, we show that even though the two- 
point mass correlation function decouples for any p, product measure is not valid for p < 1. This means that the 
correlations between masses at different sites show up only in 3 or higher order correlation functions but not at the 
2 point level. Exact calculation of the 3-point correlation function will be presented elsewhere [p3|. Interestingly 
however the single point mass distribution P{m) obtained using product measure ansatz is extremely close to the 
numerically obtained distribution. 

Interpreting rm as the distance between two hard core particles labelled i and (i + 1), it is easy to see that within 
product measure ansatz, the steady state probability of a given configuration can be written as, P{mi,m2, ■ ■ ■) ^ 
Y[i''^i~ for small gaps between neighbouring particles. This represents a gas of particles moving on a ring 
with nearest neighbour interaction P{p) log(r) for small r, where r is the separation between neighbouring particles. 
Choosing different dynamics via tuning p corresponds to changing the coupling continuously. For the asymmetric 
model, (3{p) = (1 — p)/{2 — p) for p < 1 and —1 for p = 1. This corresponds to a shift from a potential that prefers 
"bunching" of particles for p < 1 to a repulsive one at p = 1. This jump discontinuity is lost for the symmetric model 
where we have f3{p) — [2 — p)^/(8 — 5p + p^) > for all p. 

We also calculated exactly the correlation function Cj =< ttt-qTOj > for the asymmetric and the symmetric models. 
When the transport is asymmetric the correlation function factorises for j ^ 0. Unlike the asymmetric case, there 
are nontrivial correlations in the symmetric model. The connected part of the correlation function oscillates with 
distance and the amplitude of the oscillation decays exponentially with distance (see figure 5). 

A simple lattice mass model with diffusion, aggregation and chipping of single units of mass was shown to exhibit 
nonequilibrium phase transition in the steady state l3Jy,0|. In this paper we have shown that if a random fraction 
chips off instead of a single unit, the steady state no longer has a phase transition as the rates of microscopic processes 
are varied. 

There are several directions for future work. For the asymmetric continuous mass model with continuous dynamics 
(p — > limit of our model), Krug and Garcia pq ] had derived density-density correlations between particles in the 
lattice gas representation. It would be interesting to extend their calculation to general p. Another interesting 
direction would be to derive the large scale hydrodynamics for general p and extend the calculation of the tracer 
diffusion coefficient |lqj23] to general p. 

We thank D. Dhar, M. Barma and J. Krug for very useful discussions. 



APPENDIX A: NON EXACTNESS OF PRODUCT MEASURE ANSATZ FOR ASYMMETRIC 

CONTINUOUS MODEL FOR p< 1 

In this appendix, we show that the product measure or the mean field theory is not exact for asymmetric continuous 
mass model for p 7^ 1. The steps in the proof are as follows. In general the evolution equation of the n-point joint 
distribution function will involve the {n + l)-point distribution. If the product measure were to be exact, then every 
equation of this hierarchy has to be consistent with that ansatz. We show below that for the asymmetric case, the 
second equation of the hierarchy namely the one involving the 2-point and 3-point distributions is not consistent with 
product measure ansatz. 

Firstly, we recall that we can derive an expression for the single point distribution P[m) in the steady state by 
assuming product measure in the equation involving the single point and two point distributions (namely Eq. (H)). 
The Laplace transform (3(s) = pV {pV -\- q) /{I — pqV — q^) is given by Eq. (||) where (5(s) = d{sV)/ds. Next we write 
down the second equation of the hierarchy, namely the evolution equation of the joint mass distribution, P{mi, m^+i) 
of two adjacent sites i and (i + 1), 

P{mi,mi+i,t + l) = / dnii^i / dri^if{ri^i) / dm[ / drif{ri) / dm^+i / (ir,;+i/(ri+i) 

X P(mi_i,mi,TOi+i,t)(5(mi_iri_i +m^(l - n) - rrii) 

X 5{m'^ri + m'^^-^{l - r^+i) - to,;+i). (33) 

All the integrals over dm run from to 00 while the integrals over dr run from to 1. P{mi^i,mi,mi^i,t) is the 
three point joint mass distribution function and /(r) is given by Eq. (|g). If the product measure were exact, the 
joint distributions in the above equation would factorize and the resulting equation must be satisfied by the P{m) 
obtained from the first equation of the hierarchy, namely Eq. (p[ ) . 

Assuming factorization P{mi,m,2 ■ ■ ) — Yii ^("^j) ii^ Eq. (133^ , multiplying both sides by e~'"»^i~™*+i*2 g^j^^j then 
integrating over rrii and rrii+i, we get 

Q{si)Q{s2) = {q+pV{si)) {qQ{s2) +pV{s2)) (qQ{si)+pJ dr,Q{si + (53 - S2)n)] . (34) 

If product measure were to be true, this equation must be satisfied with the expression of Q{s) obtained from Eq. 
(g). If we substitute the expression for Q(s) from Eq. (|^) in the above equation, wc find after a somewhat tedious 
but straightforward algebra, that Eq. ( P4|) reduces to, 

^^^^^^"' (2-p)S^ + l-p = 0. (35) 



V{si)J ' "Visi 

If product measure is to be true then a necessary condition (but not sufficient) is that the above equation be satisfied 
for arbitrary values of si and S2- For p < 1, this is an algebraic equation for the ratio V{s2)/V{si). Since the 
coefficients do not involve si or S2, the solution for V{s2)/V{si) will be a constant independent of si and S2- Clearly 
this can not be true for arbitrary values of si and S2- Thus product ansatz is not exact for p < 1 for asymmetric 
dynamics. 

Note however that for p = 1 , Eq. ( p5| ) becomes an identity. This however is a necessary but not sufficient condition 
to prove that product measure is exact for p — 1. However it was shown p5| that for p = 1, all equations of the 
hierarchy of distribution functions are actually consistent with product measure ansatz. 

APPENDIX B: PROOF OF EXACTNESS OF PRODUCT ANSATZ FOR p^ FOR THE SYMMETRIC 

MODEL 

In this appendix we show that the mean field is exact for the p ^ limit of the symmetric continuous model. 
Consider a cluster of n consecutive sites 1,2, ... ,n. In the steady state, the joint probability distribution function 
P(mi, TO2, . . . , rrin) satisfies the equation, 

= -(2n + 2)P(mi,...,m„) 

+ / dniQ I dr I drn[P{rnQ, . . . ,rnn){S{rn'i + rnor -- nii) + S{rn'i{l ~ r) — rni)) 
Jo Jq Jo 
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/•OO pi pOO 

+ dnin+i dr (iTO'„P(mi, . . . , TO„+i)((5(rn^ + m„+ir - m„) + 5(toJj(1 - r) - m„)) 

Jo Jo Jo 

j^i "'0 "'0 Jo 

X {5{m[{\ — r) — mi)5{m[j^-^ + m[r — mi+i) + 5{m[ + m-^;^r — to,;)(5(to^^i(1 — r) — rrii+i)). (36) 

The first term is the total rate of going out of the state. The second and third terms describe the mass transfer at 
the boundary of the n-cluster while the last term accounts for mass transfer within the cluster. Let P(mi, . . . , 77i„) = 
YYi P{fni)- We multiply both sides of the equation by g"™!^!"-"™"*" and sum over mi, . . . , ?n„. The resulting terms 



in the right hand side can be simplified by using the explicit expression of V{s) from Eq. (10). Then each one of 
the terms involving the integrals reduces to 2]^"Q(si), where Q{s) = J^ P{m)e~"^'^dm as before. Thus Eq. ( p^ ) is 
indeed satisfied by the product measure ansatz for all n. Joint probability distributions for any n arbitrary sites can 
be split up into product of distributions for clusters of neighbouring sites, and then the proof can be applied for each 
of the individual clusters. 

We note that for the p ^ limit of the symmetric model, it was shown in Ref. jl^l by a different method that the 
product measure is exact for any finite system of size N. 

For symmetric model with p > 0, the product measure is not exact as was shown in the text by explicit calculation 
of two point mass correlation function. 
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FIG. 1. Asymmetric continuous mass model: A random fraction r; of each mass rrii is broken off and added to the right 
neighbour with probability p. With probability (1 — p), the broken piece rejoins the original mass. 
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FIG. 2. Mapping of the mass model to a particle model is shown. Each transfer of mass to the right corresponds to the 
particle (filled circles) jumping to the left. 
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FIG. 3. The analytical mean field answer for I3{p) is compared with the numerical result for the symmetric continuous mass 
model. While the numerical single site distribution for masses is quite different from the mean field answer (see Table II), the 
small m behaviour is predicted well by the mean field. 
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FIG. 4. The analytical mean field answer for the single site mass distribution function P{m) given by Eq. (pCJ) is compared 
to the numerical result for the asymmetric random sequential discrete mass model. The data is for system size L=20000. The 
closeness of the two curves suggest that mean field P{m) is exact even though the product measure fails. 
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FIG. 5. The correlation function C{x) =< mom^ > — < m >^ given by Eq. 
distance with the amplitude decaying exponentially to zero. 



fj is shown for /ii — 0.96. It oscillates with 



TABLE I. Numerically obtained moments of the meiss are compared with the mean field values (Eq. [l3[) for p = 0.8 in the 
asymmetric continuous model. The excellent agreement with the mean field results in this case is to be contrasted with rather 
poor agreement with mean field results in case of symmetric continuous model (see Table II) . 
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Moments 


Numerical 


Mean Field 


< m^ > 


1.7998(0) 


1.80 


<m-' > 


4.6826(0) 


4.68 


< m* > 


15.9888(3) 


15.96 


<m^ > 


67.72(1) 


67.50 



TABLE II. Comparison of numerically obtained moments of the mass with the mean field values for p = 0.8 in the symmetric 
continuous model. This clearly shows that the mean field approximation is not good for the symmetric case as compared to 
the asymmetric case. 



Moments 


Numerical 


Mean Field 


< m^ > 


2.3237(4) 


2.100 


< m^ > 


8.623(5) 


6.660 


< m* > 


44.37(7) 


28.260 


< m" > 


293.2(9) 


150.314 
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